home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / zlaein.z / zlaein
Text File  |  1996-03-14  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. ZZZZLLLLAAAAEEEEIIIINNNN((((3333FFFF))))                                                          ZZZZLLLLAAAAEEEEIIIINNNN((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      ZLAEIN - use inverse iteration to find a right or left eigenvector
  10.      corresponding to the eigenvalue W of a complex upper Hessenberg matrix H
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE ZLAEIN( RIGHTV, NOINIT, N, H, LDH, W, V, B, LDB, RWORK, EPS3,
  14.                         SMLNUM, INFO )
  15.  
  16.          LOGICAL        NOINIT, RIGHTV
  17.  
  18.          INTEGER        INFO, LDB, LDH, N
  19.  
  20.          DOUBLE         PRECISION EPS3, SMLNUM
  21.  
  22.          COMPLEX*16     W
  23.  
  24.          DOUBLE         PRECISION RWORK( * )
  25.  
  26.          COMPLEX*16     B( LDB, * ), H( LDH, * ), V( * )
  27.  
  28. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  29.      ZLAEIN uses inverse iteration to find a right or left eigenvector
  30.      corresponding to the eigenvalue W of a complex upper Hessenberg matrix H.
  31.  
  32.  
  33. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  34.      RIGHTV   (input) LOGICAL
  35.               = .TRUE. : compute right eigenvector;
  36.               = .FALSE.: compute left eigenvector.
  37.  
  38.      NOINIT   (input) LOGICAL
  39.               = .TRUE. : no initial vector supplied in V
  40.               = .FALSE.: initial vector supplied in V.
  41.  
  42.      N       (input) INTEGER
  43.              The order of the matrix H.  N >= 0.
  44.  
  45.      H       (input) COMPLEX*16 array, dimension (LDH,N)
  46.              The upper Hessenberg matrix H.
  47.  
  48.      LDH     (input) INTEGER
  49.              The leading dimension of the array H.  LDH >= max(1,N).
  50.  
  51.      W       (input) COMPLEX*16
  52.              The eigenvalue of H whose corresponding right or left eigenvector
  53.              is to be computed.
  54.  
  55.      V       (input/output) COMPLEX*16 array, dimension (N)
  56.              On entry, if NOINIT = .FALSE., V must contain a starting vector
  57.              for inverse iteration; otherwise V need not be set.  On exit, V
  58.              contains the computed eigenvector, normalized so that the
  59.              component of largest magnitude has magnitude 1; here the
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. ZZZZLLLLAAAAEEEEIIIINNNN((((3333FFFF))))                                                          ZZZZLLLLAAAAEEEEIIIINNNN((((3333FFFF))))
  71.  
  72.  
  73.  
  74.              magnitude of a complex number (x,y) is taken to be |x| + |y|.
  75.  
  76.      B       (workspace) COMPLEX*16 array, dimension (LDB,N)
  77.  
  78.      LDB     (input) INTEGER
  79.              The leading dimension of the array B.  LDB >= max(1,N).
  80.  
  81.      RWORK   (workspace) DOUBLE PRECISION array, dimension (N)
  82.  
  83.      EPS3    (input) DOUBLE PRECISION
  84.              A small machine-dependent value which is used to perturb close
  85.              eigenvalues, and to replace zero pivots.
  86.  
  87.      SMLNUM  (input) DOUBLE PRECISION
  88.              A machine-dependent value close to the underflow threshold.
  89.  
  90.      INFO    (output) INTEGER
  91.              = 0:  successful exit
  92.              = 1:  inverse iteration did not converge; V is set to the last
  93.              iterate.
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.